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Abstract 



The exact solution of the asymmetric exclusion problem and several 
of its generalizations is obtained by a matrix product ansatz. Due to the 
similarity of the master equation and the Schrodinger equation at imaginary 
times the solution of these problems reduces to the diagonalization of a one 
dimensional quantum Hamiltonian. We present initially the solution of the 
problem when an arbitrary mixture of molecules, each of then having an arbi- 
trary size {s = 0, 1, 2, . . .) in units of lattice spacing, diffuses asymmetrically 
on the lattice. The solution of the more general problem where we have the 
diffusion of particles belonging to N distinct class of particles (c = 1, . . . , N)^ 
with hierarchical order, and arbitrary sizes is also solved. Our matrix prod- 
uct ansatz asserts that the amplitudes of an arbitrary eigenfunction of the 
associated quantum Hamiltonian can be expressed by a product of matrices. 
The algebraic properties of the matrices defining the ansatz depend on the 
particular associated Hamiltonian. The absence of contradictions in the al- 



gebraic relations defining the algebra ensures the exact integrability of the 
model. In the case of particles distributed in A'^ > 2 classes, the associativity 
of the above algebra implies the Yang-Baxter relations of the exact integrable 
model. 

KEY WORDS: Asymmetric diffusion, Quantum chains, Matrix product 
ansatz, Bethe ansatz 
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I. INTRODUCTION 



The representation of interacting stochastic particle dynamics in terms of quantum spin 
systems produced interesting and fruitful interchanges among the area of equilibrium and 
noncquilibrium statistical mechanics. The connection among these areas follows from the 
similarity between the master equation describing the time-fluctuations on the nonequilib- 
rium stochastic problem and the quantum fluctuations of the equilibrium quantum spin 
chains [1]- [20]. 

Unlike the area of nonequilibrium interacting systems, where very few models are fully 
solvable, there exists a huge family of quantum chains appearing in equilibrium problems 
that are exactly integrable. The machinery that allows the exact solutions of these quantum 
chains comes from the Bethe ansatz on its several formulations (see [21]- [24] for reviews). 
The above mentioned mathematical connection between equilibrium and nonequilibrium 
revealed that some quantum chains related to interacting stochastic problems are exactly 
solvable through the Bethe ansatz. The simplest example is the problem of asymmetric 
diffusion of hard-core particles on the one dimensional lattice (see [16,17,20] for reviews). The 
time fluctuations of this last model is governed by a time evolution operator that coincides 
with the exact integrable anisotropic Heisenberg chain, or the so called, XXZ quantum 
chain, in its ferromagnetically ordered regime. A generalization of this stochastic problem 
where exact integrability is also known [25]- [27] is the case where there exists N {N = 
1,2,...) classes of particles hierarchically ordered diffusing asymmetrically on the lattice. 
The quantum chain related to this problem is known in the literature as the anisotropic 
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Sutherland model [28] or SU(3) Perk-Schultz model [29]. In [15], [18] and [19] it was shown 
that the above mentioned asymmetric exclusion problem could also be solved exactly though 
the Bethe ansatz in the cases where the particles diffusing on the lattice have hard-core 
interactions of arbitrary range, or equivalently, the particles have size s = 0, 1, 2, in units 
of lattice spacing. 

On the other hand, along the last decade it has been shown that the stationary dis- 
tribution of probability densities of some stochastic models can also be expressed in terms 
of a matrix product ansatz. This means that the ground state eigenvector of the related 
quantum chain is also given by a matrix product ansatz. According to this ansatz the com- 
ponents of the ground state wavefunction are given in terms of a product of matrices. These 
components, apart from an overall normalization constant, are fixed by the commutation 
relations of the matrices defining the matrix product ansatz. These models are in general not 
exact integrable [13] and the matrix product ansatz only gives the ground state wavefunc- 
tion of the related Hamiltonian. Despite of this limitation this ansatz produced interesting 
results in a quite variety of problems including interface growth [30], boundary induce phase 
transitions [31]- [34] , the dynamics of shocks [35] or traffic flow [36] . 

An interesting development of the matrix product ansatz that happened also in the area 
of interacting stochastic models is nowadays known as the dynamical matrix ansatz [37]. 
According to this new ansatz, whenever it is valid, the probability density of the stochastic 
system is given by a matrix product ansatz not only at the stationary state but at arbi- 
trary times. In the related quantum chain this would be equivalent to the requirement 
that not only the ground state wavefunction, but an arbitrary one, should have its compo- 
nents given by a matrix product ansatz. The dynamical matrix product ansatz was applied 
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originally on the problem of asymmetric diffusion of particles on the lattice [37,38]. More 
recently [39] , [40] it was also shown that this ansatz can also be formulated in the problem of 
asymmetric diffusion of two types of particles. The validity of the ansatz was confirmed in 
the regions where the model is know to be exact integrable through the Bethe ansatz [3,18]. 
Motivated by this fact we decided to verify if we can solve the above quantum chains directly 
though a matrix product ansatz, without considering any time dependence as happens in 
the dynamical matrix ansatz. Surprisingly we were able to rederive all the results previously 
obtained though the Bethe ansatz for the asymmetric diffusion problem with one specie of 
particles [15] or more [18,19]. Moreover our derivation turns out to be quite simple and it is 
not difficult to extend it to many other quantum Hamiltonians related or not to stochastic 
particle dynamics [41]. We are going to present in this paper these derivations and as we 
shall see, many of the results obtained in [15], [18] and [19] can now be rederived quite easily. 
The simphcity of our ansatz enabled us to extend the results of [19] to the case where each 
individual particle i belonging to any class {c — 1, . . . ,N) is distinguishable with a given 
size Si {si = 0, 1, . . .). 

The paper is organized as follows. In the next section we review the asymmetric diffusion 
problem with a single type of particles of arbitrary size and we derive the associated quantum 
chain. In section 3 we introduce the matrix product ansatz and obtain the exact solution 
of the model presented in section 2. Similarly as in section 2, in section 4, we derive the 
quantum Hamiltonian associated to the problem of asymmetric diffusion of several types of 
particles with arbitrary sizes and hierarchical order. In section 5 the exact solution of the 
general model of section 4 is obtained though an appropriate matrix product ansatz. Finally 
in section 6 we conclude our paper with some final comments and conclusions. 
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II. THE ASYMMETRIC EXCLUSION MODEL WITH PARTICLES OF 



ARBITRARY SIZES 

The standard asymmetric exclusion model is a one-dimensional stochastic model that 
describes the time fluctuations of hard-core particles diffusing asymmetrically on the lattice. 
If we denote an occupied site i on the lattice by a? — +1 and a vacant site i hy af — —1, the 
time evolution operator of the probability distribution of particles is given by the following 
asymmetric XXZ Hamiltonian 



L r 



i=l 



e+a, a, 



+ e_a+ar,, + -(l-aX+i) 



(1) 



where L is the number of lattice sizes and cr^ = ((7^±i(7^)/2 are the raising and lowering 
spin— i Pauli operators. Periodic boundary conditions are imposed and e+ and e_ (e_|_ -|-e_ = 
1) are the transition probabilities for the motions to the right and left, respectively. It is 
important to notice that this Hamiltonian, contrary to the standard XXZ quantum chain, is 
not Hermitian for e+ 7^ e_. Such property, besides producing complex value eigenvalues also 
produces massless regime, in a region where the standard XXZ is massive (gapped), whose 
mass gap vanishes as L-3/2 [10]- [12]. 

The generalization of this problem, that we consider in this section, is obtained by 
considering each distinct particle, instead of having an excluded hard-core volume of a 
single lattice size (s = 1), may now have a hard-core volume of s sites (s = 0,1,2,...). 
Equivalently, each individual particle on the lattice may have a distinct size s — 0,1,2,.... 
Particles of sizes s on the lattice are composed by s one-site monomerers and we represent 
their coordinates by giving the position of their leftmost monomer. In Fig. 1 some examples 
were shown for the configurations with n — 5 molecules and some size distributions {s} in 
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a lattice with L — 5 sites. We should notice that molecules of size s = has no excluded 
volume interaction and we can have an arbitrary number of them in a given site. However 
we should stress that although being sizeless they keep the order of the size's distribution 
on the lattice. This means that if a given particle of size s is initially between particles of 
size s' and s" it will keep this relative order in future times. 

In order to describe the occupancy of a given site i {1,2,...,N) we attach on it a site 
variable Pi taking integer values e Z). If /3j = 0, the site is vacant (or may be occupied 
by a monomer of the molecule on its leftmost neighboring site). If /9j > 0, we have on the 
site a molecule of size s — (3i and the sites j = + — 1 are empty sites. Finally 

if (5i — —n < we have, at the site i, n molecules of size zero. The allowed configurations, 
denoted by = /32, Pn} are those satisfying the hard-core constraints imposed by 
the sizes of the molecules on the periodic lattice. This means that if in a given configuration 
{13 i} we have I3j ^ and Pi ^ then we should have I — j>Piorj — l> Pj depending if 
I > j or I < j, respectively (see Fig. 1). 

The master equation for the probability distribution at a given time t, P{{P},t), can be 
written in general as 

dP{{P},t) 



at 



= -rm ^ {p'})P{{p}, t) + r{{p'} ^ {P})P{{P'}, t) (2) 



where r({/3} — > {/?'}) is the transition rate where a configuration {P} changes to {P'}. In 
the model we are considering there exists only diffusion processes. The allowed motions, 
whenever there is no hard-core constraints, are those in which a given particle diffuses to its 
right, 

A 0i+l ^ 0^ Pi+1, P>0 



A 7m^(/3 + l)i (7-l)m, /3 < 0, 7 < 0, (3) 
with transition rate er, and diffusion to the left, 

7, (7-1), (/3 + !),+!, /3<0, 7<0, (4) 

with transition rate e^. The master equation (2) can be written as a Schrodinger equation 
in Euchdean time (see [3] for general apphcations for two-body processes) 

d\P) 



dt 



-H\P), (5) 



if we interpret \P) = P({/3},t) as the associated wave function. If we represent Pi as \j3)i, 
the vectors ® \/3)2 'S' ■ • ■ ® \P)n will span the associated Hilbert space. The diffusion 
process given in (3) and (4) give us the Hamiltonian [3] 

i=l 

00 
13=1 

/3=-<x> 7=-c» 

with 

D^eR + eL, e+ = — ^ — , e_ = — ^ — , (7) 
^R + eL eR + cl 

and periodic boundary conditions. The matrices E°''^ are infinite-dimensional with a single 
nonzero element = Sa,iS/3,j {a,(3,i,j e Z). The projector V projects out the con- 

figurations > satisfying the constraint that for all Pi,Pj 7^ : (j — > Sj if j > i or 

{i — j) > Sj if i > j. The constant D in (6) fixes the time scale and for simplicity we chose 



D — 1. A simplification of our general problem happens when all the particles have the 
same size s > 0. In this case the matrices E"''^ can be replaced by the spin-| Pauli matrices 
and the Hamiltonian is given by 

where now Vs projects out the configuration where two up spins, in the cr'^-basis, are at 
distance smaller than the size s > of the particles. The simplest case s = 1 gives Ps = 1 and 
we obtain the standard asymmetric exclusion Hamiltonian (1). For the sake of comparison 
with the standard XXZ chain, normally considered in the context of magnetic systems, we 
perform for e+, e_ the following canonical transformation: 

4 (-)^^^f , ^ a^ ^ = 1, 2, . . . , L, (9) 



+ 



in (8) and obtain 



H = -\v^j:Vs + a^a^^, + A{atat^, - 1)] P„ (10) 



A 



1=1 
e+ + e_ 



Apart from the projector this Hamiltonian coincides with the gapped ferromagnetic Heisen- 
berg chain. However now, in distinction with (8), the boundary condition is not periodic 
but twisted 



a 



N+l 



''^f^at,a%^, = al (11) 



Since 7^ 1 this boundary term has the same degree of importance than the whole system 
and we have a critical behavior induced by the surface, i. e., the mass gap vanishes in 
opposition to the standard periodic ferromagnetic XXZ chain. 



III. THE EXACT SOLUTION OF THE GENERALIZED ASYMMETRIC 
EXCLUSION PROBLEM: THE MATRIX PRODUCT ANSATZ 



The exact solution of the generalized asymmetric exclusion problem of last section was 
derived in [15] within the framework of the coordinate Bethe ansatz. In this section we are 
going to rederive this solution by imposing a matrix product ansatz for the eigenfunctions 
of the Hamiltonian (6). As we shall see this derivation turns out to be more direct than the 
old one presented in [15]. 

Before considering the more general situation where any molecule may have a distinct 
size let us consider initially the simple case where all the molecules have the same size s 
(s = 0,l,...). 

Since the diffusion process conserve particles, and the lattice is periodic, the total number 
of particles n and the momentum P are good quantum numbers. Consequently the Hilbert 
space associated to (6) can be separated into block disjoints sectors labelled by the values 
of n (n = 0, 1, . . .) and P (P = ^; / = 0, 12, . . . , L - 1). 

Our ansatz asserts that any eigenfunction |^„,p) of (6) in the sector with n particles and 

momentum P, will have its components given by the matrix product 

* 

|*n,p) = E f{xi,...,Xn)\xi,...,Xn), (12) 

{xi,...,x„} 

The ket \xi, . . . , Xn) denotes the configuration where the particles are located at (xi, . . . , Xn) 
and the symbol (*) in the sum denotes the restriction to the sets satisfying the hard-core 
exclusion due to the size s of the particles, i. e., 

Xi+i > Xi + s, i — 1, . . . ,n — 1, s < Xn — Xi < L — s, (13) 
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where we have to remember that in the case where the particles have size s = we may have 
any number of particle in a given site. Differently from the standard Bethe ansatz where 
/(xi, . . . , Xn) is given by the combination of plane waves now it is given by the trace of a 
product of matrices. The matrices E and A^^"* are associated to the empty and occupied 
sites describing the configuration of the lattice. The superscript (s) is just to remember the 
size of the particle. The matrix Qp in (12) is introduced in order to ensure the momentum 
P of the eigenfunction l^'^ p). This is accomplished by imposing the commutation relation 

E^p = e-'^VtpE, A'^p = e-'■^^p^^ (14) 

since from (12) we must have for eigenfunctions of momentum P the ratio of the amplitudes: 

e--^ (m=l,2,...,L-l). (15) 



f{xi, . . . ,Xn) -imP 



f{xi + m,...,Xn + m) 

The algebraic properties of A^'^^ and E will be fixed by requiring that l^'n.p); defining the 
ansatz (12), satisfy the eigenvalue equation 

i/{.i=...=s„=.}|^'n,p) = en\'^n,p), (16) 

where i?{si=...=s„=4 is given by (8). 

Before considering the case where n is general let us consider initially the cases where 
we have only n = 1 or n = 2 particles. 

n = 1. For one particle the eigenvalue equation (16) give us 

£iTr (£;^i-U(^)£;^-^iOp) = -e+Tr (^E'^'-^A^'^ E^-'^'+^flp) 
-e_Tr (^E'^'A^'^E^-'^'-^Qp) + (e+ + e_)Tr (e^^-^A^'^E^-^^I^p) . (17) 

The cyclic property of the trace and the algebra (14) fix the values of the energies 
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£i = -(e+e-^^ + e_e^^-l), (18) 

where p — ^ (Z = 0, 1, . . . , L — 1), is the momentum of the state. 

An alternative way to solve (17) that will be easier to generalize for arbitrary values of 
n is obtained by the replacement 

- (19) 

where is a spectral parameter dependent matrix with the following commutation relation 
with the matrix E: 

EAi'^ = e^'Ai'^E. (20) 
Inserting (19) in (17) and using (20) we obtain 

£^ = £{k) = -{e+e-''' + e_e*^' - 1), (21) 

where we have used e_|_ + e_ = 1. 

Comparing (21) with (18) we fix the spectral parameter k as the momentum of the 
1-particle eigenfunction |*i,p), i. e., k — P — ^ {I — 0,1, . . . , L — 1). 

n —2. For two particles on the lattice the eigenvalue equation (16) gives for |^'2,p) 
two types of relations depending on the relative location of the particles. The amplitudes 
corresponding to the configuration \xi, X2) where X2 > Xi + s will give the relation 

-e_Tr [E''^A^'^E"^-'^A^'^E^-''^np) - e+Tr (E''^-^A^'^E"^-^A^'^E^-''^+^Vtp) 
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A possible and convenient way to solve this equation is by identifying the matrices A^^^ as 

(s) (s) 

composed by two spectral parameter dependent new matrices A^^ and A]^^, i. e., 

A^'^=j:4lE'-', (23) 

i=l 

that satisfy, as in (20) the commutation relation 

E4!-e''^Ai^E, a = 1,2). (24) 
Inserting (23) in (22) and using (24) we obtain 

S2 = e{h) + s{k2), (25) 

where e{k) is given in (21). 

The relation (14) give the commutation of these new matrices A^f,'' with Qp, i. e., 

^.^Op^e^^MOp^), (j = l,2). (26) 

Comparing the components of the configurations \xi, X2) and \xi + 'm,X2 + m), and exploring 
the cyclic invariance of the trace we obtain 

P = ki + k2. (27) 

Up to now the commutation relations of the matrices A^j^^ and A^*'' among themselves as 
well the spectral parameters, that in general may be complex, are unknown. The eigenvalue 
equation (16) when applied to the components of the configuration \xi, X2) where X2 — xi-\-s 
("matching" conditions) will give us the relation 

£2Tr {E'''-^A^'^E'-^A^'^E^-''^-'ilp) = -e+Tr (e''^-'^A^'^E'A^'^E^-'''-'^p) 
-e_Tr (£;^i-M(^)EM(^)^^-^i-*-^Op) + Tr (^^^-^^W^^-^^W^^-^i-^Op) . (28) 
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Using (23) to express the A^*) matrix in terms of the spectral parameter matrices A^^^ 
(j = 1,2), and (21) for £2, the last expression gives 

[e- - e-""^ + e+e-^(*^^+^'^] Ai'^Aif = 0. (29) 

This last relation imply that the matrices {A^. } should obey the algebra 

where 

S(ki, ki) = - -a^,, • (31) 

The complex spectral parameters {/Cj} are still free up to now. The cyclic property of the 
trace together with the algebraic relations (14), (24), (26) and (30) and the fact that any 
component should be uniquely related give us 

Tr {A^^A^.E^^-'^'^'^^p) = e-*(^-2^+2)fc.Tr (4f E^-^^+^^Jf.^f^p) 

_ g-ife,Lgi2fc,(.-i)g-iP(.-i)^(^ (4^)^(^)£;^-2^+2j^p) , (32) 

or equivalently, since P = ki + k2, 

(ik \ 
. i = l>2 (jVO- (33) 

The energy and momentum are given by inserting the solution of (33) into (25) and (27), 
respectively. 

General n. The above calculation can easily be extended to the case where n> 2. The 
eigenvalue equation (16) when applied to the components of the eigenfunction corresponding 
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to the configuration of |^n,p) where all the particles are at distances larger than the size s 
of the particles, gives a generalization of (22): 

n 
i=l 

+e_Tr (• • • _£;^i-^i-i-i^(«)_£;^i+i-^i-2^(^) . . . >l(s)£;^-^"+iOp) 

— Tr ^- • • ^^i+<-^i-i^(«) . . . ^(«) £'-'"~^"Qp^ ].. (34) 

The solution is obtained by identifying the A^^^^ matrix as a combination of n spectral 
parameter dependent {^fc*^} matrices, namely, 

^^^^ = ^:4?^'"^ (35) 

with the commutation relations with the matrices E and Qp 

= ^'^A^E, A^^^p = e^^(^-)0p4) (j = 1, . . . , n). (36) 

Inserting (35) into (34) and using the relations (36), together with the cyclic property of the 
trace we obtain 

n n 

£„ = ^£(/c,), P = E%> (37) 

for the energy and momentum of |^n,p), respctively. The eigenvalue equation (16) applied to 
the configuration where a pair of particles located at xi and Xj+i are at the closest position, 
i. e., Xj+i = + s, will give relations that coincides with (27) and (31), but now with 
j — l,2,...,n. The configurations of \^n,p) corresponding to three or more particles at 
the "matching" distances will demand that the algebra satisfied by the matrices {A^.} in 
(30) is associative. Equivalently this means that a given component, expressed in terms of a 
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product of matrices {^Ife^ } and should be uniquely related to the other components. This 



is immediate for the present problem since the structure constants S{kj,ki) of the algebra 
in (30) are constants with the property 



As we are going to see in section 4, this condition in general leads to the well known Yang- 
Baxter relations [42,21]. 

Once the algebra is defined all the components of |^n,p) can be uniquely determined 
only if this algebra has a well defined trace, whose cyclic property will fix the n complex 
spectral parameters {%}. An analog procedure as in (32) give us the constraints 



This equation coincides with the Bethe-ansatz equations derived in [15] through the coor- 
dinate Bethe ansatz method. Moreover, an arbitrary component f{xi, . . . , x„) of the wave 
function |^n,p) given in (12), can be written as 

fi^li ■ ■ ■ 1 ^n) — 



S{kj,ki)S{ki,kj) = 1. 



(38) 




(39) 



n 



n 



n 




(40) 



Using the commutation relation (24) and the fact that (^i^ j = (j = 1, . . . , n) 



we can 



rewrite this last expression as 



fipl-i • • • -I ^n) 




(41) 



Let us define the new matrices 
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(42) 



It is simple to verify, from (36), that they satisfy 



(43) 



where 




Finally, in terms of these new matrices, and exploring the fact that (^i^^) = we can write 



f(x„ . . . , x„) = E e^[*^..(-i-i)+^..(-.-i)+-+^.n(^.-i)]rr . . . AiiE^np) , (45) 



Pl,--;Pn 



where the sum is over the permutations (pi,P2, • • • ,Pn) of non repeated integers (1, 2, . . . , n). 
The result (45) show us that the amplitudes derived using the present matrix product ansatz 
is given by a combination of planes waves with complex wave number {kj}, and reproduces 
the results obtained previously [15] through the standard coordinate Bethe ansatz. 

Let us return to the general case where we have n molecules with arbitrary sizes 
{si, S2, ■ ■ ■ , Sn} and whose related Hamiltonian is given by (6). In this general case each 
particle is conserved separately, and since in the diffusion processes no interchange of par- 
ticles are allowed, also the order {si, S2, ■ ■ ■ , s„} where the particles appear is a constant 
of motion, up to cyclic permutations. The eigenfunction corresponding to a given order 
{si, S2, ■ ■ ■ , Sn} and momentum P can be written as 



where ' • (xi, . . . , Xn) is the component of a configuration where the particles of sizes 
Sci, ■ ■ ■ , Sc„ are located at positions xi, . . . ,Xn respectively. The summation {c} extends over 




(46) 



{c} {x} 
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all cyclic permutations {ci, . . . , c„} of integers {1, . . . ,n}, and the summation {x} extends, 
for a given distribution {sci, • • • , Sc^} of molecules, to increasing integers satisfying 

Xj+i — Xi '> Sf,^-, i — 1, . . . ,Tl — 1 

Sci < Xn- Xi < N - Sc„. (47) 

In order to formulate our matrix product ansatz we associate to the sites occupied by 
the particles of size sj {j = 1, . . . ,n) a matrix A^*^) and to the remaining L — n sites we 
associate, as before, the matrix E. Our ansatz asserts, in a generalization of (12), that the 
amplitudes of the eigenfunctions (46) are given by 

where in order to ensure the momentum P of the eigenstate the matrices {^'f*^} should 
satisfy 

Eflp = e-'^QpE, A^'^Qp = e-'^'flpA^'l (49) 

Let us consider initially our ansatz (46)- (49) for n = 1 and n — 2 molecules. 

n=l For one particle on the chain we have the same ansatz as (12) and the energy given 
by (18). 

n=2 In this case if both particles have the same size si — S2 — s we have the same 
situation considered previously in (22)-(33). The eigenfunctions |^I/{s,s},p) will be given by 
(12) and the energy by (24) with ki fixed by (33). 

If the particles are distinct the situation is new. The eigenvalue equation when applied to 
the configurations where the two particles of sizes Si and S2 are located at Xi and X2 > Xi+Si, 
respectively, will give, mutatis mutandis, an expression similar to (32). The corresponding 
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situation are obtained by introducing the generalization of the spectral parameter matrices 
defined on (23), i. e., 

^(s,)^^^(^.)^2-.,^ (J = l,2), (50) 
1=1 

that satisfy, as (24), the commutation relations 

^47^ = e^'^'Al;^^^, (j,/ = l,2). (51) 

The energy S{si,s2} ^nd momentum P are related to the spectral parameters by (25) and 
(27), respectively, and from (49) and (50), we have the commutation relations, generalizing 
(26) 

4f Qp = e^^(^-^)Qp4f , (j,l^l,...,n), (52) 

with n — 2. 

If the two particles are at the closest distance X2 — Xi + Si ("matching" condition) the 
expression (28) should be replaced by 

e{s,,s2}Tr [E''^-^A^'^^E'^-^A^'^^E^-''^-'^np) = -e+Tr [E''^-^A^'^^E'^A^'^^E^-'''-'^np) 
-e_Tr [e""^-^ A^'^^ E'^ A^'^^ E^-'^^-'^-^VLp) + Tr (e""^-^ A^'^^ E'^'^ A^'^^ E^'^'^-'^VLp) . 

Inserting the definition (50), the expression (25) for £{si,s2} and using the algebraic relations 
(51) and (52) we obtain the commutation relations for the matrices {^^*'^} 

4:^4:^-S{kM4:}4;\ {j^m;l,r^ 1,2), AifAif^O, (53) 

where S{kj,km) is given by the same expression as (31). It is interesting to notice that 
the structure constants S{kj,km) of the algebra in (53) are independent of the superscript 
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of the matrices A^kf, and consequently the algebra among the {^i^'^} is the same as that 
of (30) with respect to the interchange of spectral parameters. However the superscript of 
these matrices can not be neglected since in the commutations relations (53) they are not 
interchanged and also their commutation with the Jlp matrix is size dependent (see (52)). 

The spectral parameters ki and k2 are fixed by the cyclic property of the trace, and we 
have, mutatis mutandis, a similar expression as (32). Using the algebraic relations (51), (52) 
and (53) we obtain 

Tr [Ai'^'UifE''-''-''+^np) = 

^-ikjL^ik,{s,+s2-2)^-iP{s,-i)g^j^,^ j^^-^rj^ (^AifAifE^-'^-'^+^np) . (54) 

However differently from (32) the traces on the left and in the right hand side of the equation 
are not the same if si 7^ S2. If we repeat once more the commutations that lead to (54) we 
obtain the same trace in both sides and consequently 

Since P = ki + k2, this last expression is equivalent to 

e''^^ = e'^- ( 1 S{kj,ki), m = 0,l; j ^1^1,2; ^ S2, (56) 



and 



_ S1 + S2 , . 

s = — ^ — (57) 



is the average size of the two molecules. The expression (56) generalizes the expression (33) 
obtained for particles of equal sizes. We note however that since m = 0, 1 in (56) we have two 
times more solutions than the corresponding one (33) for particles of equal sizes. This indeed 
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should be the case since particles of distinct sizes are distinguishable, and consequently the 
number of eigenfunction is doubled when compared with the indistinguishable case (si = 52)- 
General n. In this case we have a general distribution of particles with sizes 
{si, S2, . . . , Sn} and the corresponding eigenfunctions are given by (46) and (48). The eigen- 
value equation when applied to the components |a;i, . . . , Xn) where all the particles are not 
at the closest distance, gives an equation similar to (34), mutatis mutandis, whose solution 
is given by the generalization of (50), (52) and (51): 

n 

A^'^^ = E4?^^'"^^ Af^P = e^^(i-^)np4f , = e^'^Ai^E, {j = 1, 2), (58) 

1=1 

producing the energy and momentum given by 

n n 

s^ = Y: (e+e-^'^^ + e_e^^^- - l) , P = k„ (59) 
3=1 j=i 

respectively The eigenvalue equation applied to the components where a pair of particles 

{xi, Xi+i) are located at the closest distance, Xi+i — Xi + Si will give a generalization of (53) 

^;^4:^ = s{k„k)Ai:;^A^;\ u^i), 4^4-) = o, {j,i,t,u = i,...,n). (eo) 

The cyclic property of the trace in (48) will give, by using (58) and (60) a generalization of 
(54), namely, for each kj, 

(n s{k,, Tr (4r^4f • • • Ai:r'^:i ■ ■ ■ ^r^E^-^t^^^-^^np) . (6i) 

Similarly as happened in (54) the traces in both sides of the last equation are not the same 
because {si, S2, ■ ■ ■ , Sn} {^n, Si, ■ ■ ■ , Sn-i}- But we can redue the above commutations by 
r times until we reach the same distribution of sizes, where r is the minimum number of 

21 



cyclic rotations of {si, S2, ■ ■ ■ , s„} where the configuration repeats the initial one. In this 
case we obtain 



II, 



1=1 



(62) 



Since P — YJi=i kj, we can rewrite this last expression as 

n / ikj \ 

e"''''-e''^^l[S{kj,k){-j-] , m = 0,l,...,r-l, j,l = l,...,n, 



(63) 



1=1 



where as in (57) 



n 



(64) 



is the average size of the molecules in the distribution {si, S2, ■ ■ ■ , Sn}- This equation gives 
us a number of solutions of order r times larger than the corresponding number in the case 
where all the particles have the same size. 

The equation (63) that fix the spectral parameters of the matrices coincides with the 
Bethe - ansat2; equations derived in [15]. Similarly as we did in (40) - (45) we can show that 
indeed the eigenfunctions we obtained by using our matrix product ansatz coincides with 
the ones derived in the framework of the coordinate Bethe ansatz. 



IV. THE ASYMMETRIC DIFFUSION MODEL WITH N CLASSES OF 
PARTICLES WITH HIERARCHICAL ORDER 

The extension of the simple exclusion problem to the case where we have N distinct 
classes of particles {c— 1,2, . . . , N) diffusing asymmetrically is not exact integrable in gen- 
eral. However the integrability of the problem can be preserved if the diffusive transitions 
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of the several species happen in an hierarchical order. This problem was considered orig- 
inally in the case = 2 as a model to describe shocks [25]- [27] in nonequilibrium. The 
stationary properties of the N — 2 [31] and N — 3 [43] models can also be studied through 
a matrix product ansatz. In [19] a generalization of this problem was considered in which 
the particles in each of the N classes {c— 1,2,. . . ,N) may have distinct sizes (si, . . . , s^), 
respectively. The solution of this generalized problem was obtained through the coordinate 
Bethe ansatz [19]. In the next section we are going to show that the solution of this problem, 
similarly as we did in the last section, can also be obtained through an appropriate matrix 
product ansatz. 

In this generalized problem we consider the particles in each class c as composed by Sc 
monomers, thus occupying Sc sites on the lattice {c— 1, 2, . . .). We consider as the position 
of the molecule the coordinate of its leftmost monomer. The excluded volume of a particle in 
class c is given by its size Sc {c — 1,2, N) in units of lattice spacing. The configurations 
of the molecules in the lattice is described by defining at each lattice site i a variable (3i 
{i = 1,2, . . . , L), taking the values /3j = 0,1, . . . , N. The values /5j = 1, . . . , represent 
sites occupied by molecules of class c = 1, 2, . . . , A^, respectively. The sites attached with 
the value /3 = are the vacant sites or those excluded due to the size of the particles. As 
an example {/?} = {1, 0, 2, 0, 2, 0} may represent the configuration where in a L = 6 sites we 
have a particle of class 1 and size si = 2 at the site 1 and two particles of class 2 and size 
S2 = 1 located at the sites 3 and 5. The allowed configurations are given, in general, by the 
set {Pi} {i — 1,2, . . . , L), where for each pair Pj) ^ with j > i we have j — i > sp^. 
The allowed stochastic process in the problem are just given by the exchange of particles 
or the asymmetric diffusion if the constraint due to the size of particles is satisfied. The 
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possible motions of a given molecule are diffusion to the right 

A 0i+i ^ 0i A+i (rate Fr), (65) 

diffusion to the left 

0i A+i ^ A 0i+i (rate T^,), (66) 
and interchange of particles 

^A'A+v (/^</^') (rate r^), 

(/3>/3') (rate r^), (67) 

with P' — 1,2, . . . , N. As we see from (67), particles of a class c interchange positions with 
those in classes c' > c with the same rate as they interchange positions with the vacant sites 
(diffusion). However the net effect of these motions is distinct from the diffusion processes, 
since by interchanging positions, distinctly from the diffusion process, the particles moves 
by Sfj lattice units, accelerating its diffusion if Sc' > 1. The identification of the master 
equation as an Schrodinger equation as in (5) will give us the Hamiltonian [19] 

H = Dj2H, 

N 

a=l 

N N 

+ E E ^«,,(^'"^.4X„ - ErE^is,E^fsJ}r, (68) 

a=l 13=1 

with 

D^Vu + Vl, e+ = -^^, e. = —^ (e+ + e_ = l), (69) 

l R-\- l L i- R-r i- L 
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' e+ a < (3 

= I a = P (70) 
, e_ a> (5 

and periodic boundary conditions. The matrices E""'^ are (A^ + 1) x (A^ + l)-dimensional ma- 
trices with a single nonzero element {E°''^)ij — SaiSpj (a, (3,1, j — . . . ,N). The projector V 
in (68) projects out the configurations associated to the vectors representing molecules 
at forbidden positions due to their finite size. Mathematically this condition means that for 
all i,j with I3j 7^ we should have \i — j\ > sp^ {j > i). The constant D in (68) fixes the 
time scale in the problem and we chose D = 1. The Hamiltonian (68) corresponding to the 
particular case where all the molecules have unit size is related to the spin-^ SU(A'" + 1) 
anisotropic Sutherland chain [28,44] or SU(A'" + 1) Perk-Schultz model [29] with twisted 
boundary conditions [3]. 

At the end of the next section we are going to present the solution of a even further 
generalized model whose solution were not derived in [19]. The solution of this model is 
quite complicated through the standard coordinate Bethe ansatz. As we shall see, nowever 
its derivation is not difficult through our matrix product ansatz. In this generalization 
instead of having all the particles in a given class c with fixed size Sc, each individual 
particle may have an arbitrary size. In this case the configurations on the lattice are given 
by {/?} = {Pi, P2, ■ ■ ■ , Pl} where /3j = (c, s) mean that the lattice site i {i = 1,2, ...,L) 
is occupied by a particle of class c (c = 1,2, ... ,N) having size s {s — 1, 2, . . .). The 
Hamiltonian related to this stochastic problem is given by a generalization of (68), namely 



L 

J 
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+ E E 



(71) 



0={c,s)0>={c>,s') 



with e+, e_ and ecc' (c, c' = 1, 2, . . . , A'") given as in (69) and (70). 

V. A MATRIX PRODUCT ANSATZ FOR THE GENERALIZED DIFFUSION 
PROBLEM WITH N CLASSES OF PARTICLES WITH HIERARCHICAL ORDER 

The exact solution of the asymmetric diffusion problem with N classes of particles, whose 
related Hamiltonian is given by (69) was obtained in [19] through the coordinate Bethe 
ansatz. In this section we are going to reobtain this solution by an appropriate matrix 
product ansatz. Moreover our solution enables the extension to the more general problem 
discussed in the last section and whose Hamiltonian was introduced in (71). Let us consider 
initially the simple case where all the particles in a given class (c = 1, . . . , A^) have a fixed 
size {sc — 1,2, . . .). Due to the conservation of particles in the diffusion and interchange 
processes the total number of particles in each class is conserved separately and we can 
split the associated Hilbert space into block disjoint eigensectors labeled by the numbers 
ni, 712, . . . , nN {ni — 0, 1, . . .) of particles on the classes i {i — 1,2, . . . , N). We want to 
obtain the eigenfunctions \ni, . . . ,11,^) of the eigenvalue equation 



H\ni, . . . , njv) = Snlni, . . . , njv). 



(72) 



where 



* 




(73) 



{c} {x} 
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The ket \xi, ci; . . . ; Xn, Cn) means the configuration where particles of class q (cj = 1, . . . , N) 
is located at position Xi {xi = 1, . . . ,L) and the total number of particles is n = ni + 
• • • + un- The summation {c} = {ci, . . . , c„} extends over all the permutations of n integer 
numbers {1,2,..., N} in which rii terms have the value i {i — 1, . . . , N), while the summation 
{x} = {xi, . . . , Xn} runs, for each permutation {c}, into the set of the nondecreasing integers 
satisfying 

Xj+i Xi Sci, i — 1, . . . ,n — 1 

Sci<Xn-Xi<L-Sc^. (74) 

The matrix product ansatz we propose asserts that an arbitrary eigenfunction 
\ni, . . . jTijsf) with momentum P, will have the amplitudes in (73) given in terms of traces of 
the matrix product 

f{xi, a;...; Xn, Cn) = Tr . . . ^x„-a:„_i-ly(c„)^L-x„^^j _ ^75-) 

The matrices Y^^'> (c = 1, . . . , A^), E and flp will obey algebraic relations that ensure the 
validity of the eigenvalue equation (72). The momentum P of the state, analogously as in 
Sec. 3 is fixed by requiring the relation 

Eflp^e-'^flpE, Y^^'>flp ^ e-'^npY^^\ c^l,...,N. (76) 

Let us consider the simplest cases of n = 1 and n — 2 particles before consider the case 
where n is general. 

n=l In this case the problem is the same as that of Sec. 3 and we obtain the energies 
given by (18). 

n=2 For two particles of classes Ci and C2 (ci, C2 = 1, . . . , A?") on the lattice we have two 
distinct types of relations depending if the amplitudes are related or not to the configurations 
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where two particles are at the closest distance X2 — xi + Sci ■ The eigenvalue equation when 
applied to the components where the particles of class Ci and C2 are at positions {x2,Xi), 
such that X2 > xi + Sci, give us the relation 



S2Tr 



^2:i-l-^(ci)^2;2-a;i-ly(c2)^L-x-2Q 



e_Tr 



e_Tr 



^xi-lyid) ^X2-xiy{c2) £iL-X2-1q^ _|_ ^jxi-lyici) ^X2-xi-1y{c2) ^L-X2Q^ _ ^YY^ 



A solution of this equation is obtained by identifying the matrices Y^'^^ as composed by two 
spectral parameter dependent new matrices Y^f and Yj^^\ i. e. , 



y(c) ^ y (c) ^2-Sc 



(78) 



i=l 



that, as in (24), satisfy the commutation relation 



3 3 



(79) 



In terms of the unknown spectral parameters kj (j = 1, 2) the energy and momentum are 
given by 



£2 = £(^1) + £(^2), P^ki + k2, 



(80) 



where e{k) — — (e+e ^'^ + e^e^^ — 1)- As a consequence of (76) and (78) we also have 



Yl^np^^^'^-^'^npYl^ (j = l,2; c=l,...,iV). 



(81) 



The eigenvalue equation (72) when applied to the components of (73) where the two particles 
are at the closest distance, i. e., X2 — xi -\- give us 
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^a;i-ly(ci)^Sci-ly(c2)^L-a;i-SciQ 



-e+Tr 



^a;i-2y(ci) ^Sciy(c2) ^L-xi-ScjQ 



- e_Tr 



— 6 TV 



^Xl - 1 y (C2 ) ^SC2 - 1 y (ci ) £jL-Xl -SC2Q 



+2(1 + eeiC2)Tr 



^xi-ly(ci)^Sci-ly(c2)^i'-a;i-Scij^ 



(82) 



Substituting (78) and (80) in this last expression and using (79) we obtain 



-i{ki+km) 



l.m 



+ e_) + e-^*^'(l - eei,c2)] Tr [£;^iy^('^i)y^(^2)^L-xi-.ci-5c2+2Q^ 



(83) 



This last equation is satisfied by imposing the following commutation relations among the 

operators {Y^''^} 



~Cl,C2j 



(84) 



where 



-(e+ + e_e^('='+'='")). 



(85) 



It is interesting to consider separately the cases where the two particles belong to the same 
class ci = C2 from the case where ci ^ C2- li ci — C2 — c {c — 1, . . . , N), since ec,c — and 
1^i,m + e*^" 7^ 'Dm,i + e*^' for / 7^ m we obtain from (84) and (85) 



(86) 



where 



(87) 



and (Z, m = 1, 2; c = 1, . . . , A?"). The relation (84) in the cases where ci 7^ C2 give us the 
equations in matrix form 
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E 

Lm,=l 



ry(ci)y(c2) 

y{c2)y(ci) 



0. 



(88) 



Similarly as in [19] the above equation can be rearranged straightforwardly by giving us the 
algebraic relations 



T^(ci)-^(C2) 



N 



r' r' —1 



(89) 



where {l,m — 1,2), ci, C2 = 1, . . . , A?" and the "structure constants" of the algebra are the 
components of an ^'-matrix whose non-zero components are given by (87) and 



SZ':cl{ku k2) = [1 - ec„c2*(^i, k2)]SZlll{k,, k^) (ci, C2 = 1, . . . , TV) 
SllfAki. k2) = e,,,,Mki, k2)S^lf,lik,, k^) (ci, C2 = 1, . . . , iV; c, ^ C2), 



(90) 



where 



(91) 



The complex parameters {ki, /C2), that are free up to now, are going to be fixed by the cyclic 
property of the trace in (75) and the algebraic relations (76), (78), (79) and (81). 

Instead of solving for the spectral parameters for the n = 2 let us consider the case of 
general n. 



General n. In this case the eigenvalue equation (72), when applied to the components of 
the eigenfunction corresponding to the configuration where all the particles are at distances 
larger than the closest distance, give us a generalization of (77) that is promptly solved by 
identifying, as in (78), the matrices Y^'^^ as combinations of n spectral parameter matrices, 
i. e., 
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Y^''^ =Y.Y^^E^-'% (92) 

i=l 

satisfying the following algebraic relations with the matrices E, 

EY^ = e^^.Y^E (j = l,...,n; c = 1, 2, . . . , TV), (93) 

and from (76) 

r^f = e'^(^--)Oprif £; (j = l,...,n; c = 1, 2, . . . , iV). (94) 



The energy and momentum in terms of the spectral parameter {kj} are given by the gener- 
alizations of (80), namely, 

n n 

Sn^^eik,), P = Y.kr (95) 

The components of the eigenfunctions corresponding to the configurations where a pair of 
particles of classes ci and C2 are located at the closest positions xi and Xj+i — Xi-\- , will 
give relations that reproduces (89) - (91). 

Since in the general case we have the product of n operators {Yj^^^}, the algebraic relations 
(76), (79), (81) and (89) should provide a unique relation among these products. For example 
the product • • • Y^fY^^^Y^]^ • • • can be related to the product • • • Y^'^^Yj^^^Y^f • • • by two 
distinct ways. Either by performing the commutations a^j /3aj jSja j^a or by 
q;/37 — > a'yP — > 7q;/3 — > 7/3q;. Consequently we should have 

N ^ 

X! (^1' (^1' ^3)'S'^','^"(^2, ^3) = 

7,7',7"=1 

x: s^r^:ik2,ks)s^;j:{k,,ks)sj%{k,,k2), (96) 

7,7',7"=1 

for a, a', a", /3, (3" — 1, . . . ,N . This last relation is just the Yang-Baxter relations [42,21] 
of the S'-matrix defined in (87) and (90). Actually the condition (96) is enough to ensure 
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that any matrix product is uniquely related and it implies the associativity of the algebra 
of the operators {Y^^^}. We can verify that the Yang-Baxter relation (96), with S given 
by (87) and (90) is satisfied by an arbitrary number of distinct species of particles A^. It is 
interesting to remark that in our solution, on comparison to that presented in [19] through 
the coordinate Bethe ansatz, has the advantage that the derived S'-matrix do not depend 
on the size of the particles and the associativity condition or Yang-Baxter relation (96) is 
easier to be verified since it is the same, independently of the particle's sizes. 

The spectral parameters {kj} are fixed by the cyclic property of the trace in (75). For 
each spectral parameter kj {j — l,...,n) the commutations relation (76), (78), (79) and 
(81) applied j times, enable us to move the operator Y^'^^'' to the left 



Tr 



y-Jci) _ _ _ y{cj.i)y{cj) _ _ _ y(cn) -^L-J^^^i^jSj+n^ 



7-1 



E E ■■■ E {e^^^^"-^'55;|(fc,,A^)...5S::::F(V2,A^)^r;:^^ 



Tr 



kj — 1 /cj-|-i kji kj P 



(97) 



Moving the operator . ^ for more n — j times to the left and using the identity 



N 

E 

c'',c'',,=l 



1 ' 1 



(98) 



we can write 



Tr 



= e 



E < ci,...,c„|r|c'i,...,4 > T: 



r' r' 
^1 v"5*-ri 



ki 



kn 



(99) 



where 



< {c}\r\{c'} >^ E {s:y^ik,,kj)--.s':f^ik,,kj)---s^^{k^,kjy^^^^^ (loo) 
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We identify in (100) T{kj;{ki}) as the N'^ x A'""-dimensional transfer matrix of an inho- 
mogeneous vertex model (inhomogeneities {h}) with Bohzmann weights given by (87) and 
(90). The model is defined on a cylinder of perimeter n with a seam along its axis producing 
the twisted boundary condition 

S2^^{K, k) = 52:;|(^„, k)e^^^^'-'\ (101) 

where as always, P is the momentum of the eigenstate. The relation (100) give us the 
conditions for the spectral parameters 

e-^'=^-(^+"-^>=i"^^*) = A(%,W) i = l,...,n, (102) 

where A(/cj, {ki}) are the eigenvalues of the transfer matrix (100). The problem of fixing the 
spectral parameters {kj} reduces to the evaluation of the eigenvalues of the transfer matrix 
(100). This can be done straightforwardly through the coordinate Bethe ansatz as in [19]. 
Extracting these eigenvalues from [19] we obtain equations that coincide with the Bethe 
ansafo equations for this general model (see Eqs. (71)-(73) in [19]). 

Let us consider the more general Hamiltonian given in (71), where now each molecule 
have an arbitrary size, independently of the sizes of the other molecules belonging to its 
class. The solution of this problem was not derived through the coordinate Bethe ansatz 
since is not simple in that formulation. The Hamiltonian (71) is composed by block disjoint 
eigenvectors labeled by {ci, Si; . . . ; c„, s„} (cj = 1, . . . , A?"; Sj = 1, 2, . . . ; j = 1, . . . , n) that 
specifies the classes and sizes of each individual particle. An arbitrary eigenfunction of (71) 
is given by a generalization of (73), namely 

|ci, Si, . . . , Cji, Sji) — ^ ^ ^ ^ fip^l-i Ci, Si, . . . , Xyi^ Cm Sfi^ l^li Ci, Si, ... , Xm C^, Sj^) , (103) 
{c,s} {x}* 
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where ci, si; . . . ; Xn, c^, s„) denotes the configuration where the particle located at Xi 
{i = 1, . . . ,L) belong to the class q (q = 1, . . . , N) and has size Sj (sj = 1, 2, . . .). The 
summations {c,s} — {cp^, Sp^; . . . ; Cp„, extends over all the permutations of particles. 
The summation {x}* — {xi, . . . , Xn}* runs, for each permutation {c, s}, in the set of non- 
decreasing integers satisfying 

si < a;„ — a;i < L — (104) 

Our matrix product ansatz asserts that the eigenfunctions with a given momentum P, have 
amplitudes 

f(^Xi, Ci, <Si, . . . , Xyi, Cji, <S^) 

"Jlj. ^a;i — l-^(ci, 81)^X2— a;i — l-^(c2,S2) . . . ^a:n— a^n-l — l"J^(Cn,Sn)^i'— a^nQ^ 

where the matrices E are associated to the empty sites and F^'^j'^j) to the sites occupied by 
particles of class Cj and having a size Sj {j — 1, . . . ,n). The matrices Qp, as before, fix the 
momentum P of the eigenstates 

Eflp = e-'^QpE, r^^'^^Qp = e-'^npY^^''\ (105) 

The solution of this general problem follows, mutatis mutandis, the derivation we did in 
(77)-(102). The energy and momentum of the eigenstate is given by (95) where the spectral 
parameters are introduced by the generalization of (92): 



y(c,.) _^y^c,s)^2-.^ (106) 
i=l 

where the spectral parameter matrices Y^'^'^^ satisfy the algebra 
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EYtf^e^^YtfE = c = 1, 2, . . . , TV; s = l,2,...), (107) 

and from (105) 

Ytfnp^e^^^'-^^npYtf (j = l,...,n; c = 1, 2, . . . , TV; s = l,2,...). (108) 
The algebraic relations among the {^fc^^'*''} are given by the generalization of (89), namely, 

N 

yt'Mr^^ = 0, (109) 

where (/, m = 1, . . . , n), ci, C2 = 1, . . . , and S'^J^'^^{ki, km) are the same ^-matrix defined in 
(87), (90) and (91). It is interesting to observe that the condition of existence of a single re- 
lation among the words • • • n^"''' V^^'^^ V^'^'^^^ • • • and • • • ^^'^'^^if '''^^4"'''^ " " " reproduces 
the Yang-Baxter relation (96), as before. 

The spectral parameters {%}, as before, are fixed by the cyclic property of the trace 
and the algebraic relations (105)- (109). Using these relations we can move the operator 
Y^"''''^ to the left as in (97)- (100), 



Tr 



^ < Ci,...,Cn|'f|c'i,...,C^ > TV 



r' r' 

where now 



(110) 



<{c}\f\{d}>= {^:j;|(fci,fc,)---^S:;^(fcn,A:,)}, (HI) 

r" ••• r" 

'-'1 5 ^^n 

is different from (100) since now it corresponds to a transfer matrix of a vertex model in 
a cylinder of perimeter n with no seam (periodic boundary condition). If we iterate n — 1 
times the procedure used in obtaining (110) we obtain: 
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^ <ci,...,cjr-|c;,...,<>iv 

r' r' 



(112) 



Consequently the spectral parameters {kj} should satisfy 

e-^'=^(^+"-Sr=i^*) = e^^V^(^"-^)A(%,{A;J) j = l,...,n; r = 0, 1, . . . , n - 1, (113) 

where A(A;j, {/c;}) is an eigenvalue of the transfer matrix T given in (111), and s = Ej=i ^ 
is the average size of the particles. The eigenvalues A.{kj, {h}) can be obtained from the 
diagonalization of (111) through the coordinate Bethe ansatz and they are given by Eqs. 
(67) and (70) with = 1 of [19]. We finally have the conditions that fix the spectral 
parameters of this general problem 



X n / , ^ J = l,2,...,n, (114) 

where the auxiliary complex parameters {kf' , I — 0, . . . , N — 1; j — 1, . . . ,mi} are fixed by 
the equations 

/3=i e+ + e_e " — e s=i e+ + e_e''''^« — e'*« 

X n ^ Z = 0,l,...,7V-2; a; = l,...,m,, (115) 

where rij {j — 1,...,N), as before, is the number of particles on class j and rrii — 
J2f=i'n'j-, I — 0,...,N; mo = n,'mN = 0, and = kj. The energies and momentum 
are given in terms of {kj} by (95). We can see from (114) and (115), since for each value 
of r (r = 0, 1, . . . , n — 1) we have distinct solutions, the number of solutions is higher than 
that obtained previously. This should be expected since the particles now are completely 
distinguishable. 
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VI. CONCLUSIONS AND GENERALIZATIONS 



We have shown that all the exact results derived for the asymmetric exclusion problem 
and generalizations, through the coordinate Bethe ansatz can also be obtained in an elegant 
and unified view by an appropriate matrix product ansatz. According to this ansatz the 
amplitudes of the eigenfunctions of the associated Hamiltonian are given by traces of a 
product of matrices. The algebraic properties of the matrices appearing in the ansatz are 
fixed by the eigenvalue equation of the Hamiltonian. The existence of a well defined ratio 
among the several amplitudes of any eigenfunction imply the associativity of the algebra 
ruling the matrices defining the ansatz. In the case where we have more than a single 
kind of particles the condition of associativity of the algebra (see (96)) coincides with the 
Yang-Baxter relations [42,21]. Once the algebraic relations of the matrices are fixed the 
eigenfunctions we obtain coincides with those obtained through the coordinate Bethe ansatz. 
As an example see (45) for the case of diffusion of one kind of particle with a fixed size s. 

Differently from the Bethe annate solutions presented in [15], [18] and [19] the matrix 
product ansatz we formulate allow us to treat in a unified way the hard-core exclusion effects 
produced by the size of the particles. This virtue allowed a simple derivation of the quite 
complicated problem (see (114) and (115) in section 5) where we have N types of particles 
hierarchically ordered, but each particle being distinguishable and with a given specified 
size. The corresponding calculation through the coordinate Bethe ansatz is rater difficult. 

The extension of the solution presented in section 5 for the cases where the molecules 
are allowed to have a zero size is immediate and follows the same reasoning of sections 2 
and 3. In the case of a single specie of molecule we can also extend our models allowing 
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the molecules to have negative size (s = —1, —2, . . .) as in [45]. In this case since we do not 
have the interchange process, the particles have a well defined order on the lattice, apart of 
cyclic rotations, i. e., {xi < X2 < ■ ■ ■ < A particle i with negative size s allows a partial 
break of this ordering, i. e., {xi < X2 < ■ ■ ■ < Xi-i — s < Xi < Xj+i < • • • < Xn)- 

We may also extend the matrix product ansatz presented in this paper to the cases 
where the lattice size has open ends [41]. In those cases instead of the trace operation 
defining the amplitudes of the eigenfunction we have a single undefined matrix product, 
that can be fixed by a normalization of the corresponding eigenfunction. 

The success of our matrix product ansatz can also be tested [41] on an enormous 
variety of known exactly integrable models, irrespective if the Hamiltonian is related or 
not to nonequilibrium stochastic models. We have shown that our matrix product ansatz 
can provide the exact solutions of the XXZ chain with arbitrary exclusion effects [46] , the 
Fateev-Zamolodchikov model [47], the Izergin-Korepin model [48], the t-J model [49], the 
Hubbard model [50] as well the generalized integrable models presented in [51] and [52]. In 
conclusion our results suggest the conjecture that all exact integrable model may have its 
eigenfunctions given by an appropriate matrix product ansatz. 
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Figure Captions 



Figure 1 - Example of configurations of molecules with distinct sizes s in a lattice of size 
L = 5. 
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